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Abstract 

In this paper, we consider two-dimensional M = (4, 4) supersymmetric Yang-Mills 
(SYM) theory and deform it by a mass parameter M with keeping all supercharges. We 
further add another mass parameter m in a manner to respect two of the eight super- 
charges and put the deformed theory on a two-dimensional square lattice, on which the 
two supercharges are exactly preserved. The flat directions of scalar fields are stabilized 
due to the mass deformations, which gives discrete minima representing fuzzy spheres. We 
show in the perturbation theory that the lattice continuum limit can be taken without any 
fine tuning. Around the trivial minimum, this lattice theory serves as a non-perturbative 
definition of two-dimensional M = (4, 4) SYM theory. We also discuss that the same 
lattice theory realizes four-dimensional M = 2 U(k) SYM on R 2 x (Fuzzy R 2 ) around the 
minimum of fc-coincident fuzzy spheres. 
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1 Introduction 



Supersymmetric gauge theories play very important roles in theoretical particle physics. 
They are not only promising candidates of physics beyond the standard model, which are 
one of the targets in Large Hadron Collider (LHC) experiments, but also provide crucial 
insights into non-perturbative aspects of superstring/M theory [U |2j |3j H] . Although they 
are analytically more controllable than non-supersymmetric theories, there have arisen 
many intriguing features conjectured from various duality arguments which cannot be 
addressed by current analytic techniques. Therefore it is important to find promising 
numerical frameworks which enable us to examine them and to obtain new insights into 
non-perturbative dynamics. However, it is not a straightforward task in lattice field theory 
because of the notorious difficulties of lattice supersymmetry. So far, for one- and two- 
dimensional theories [5] and Af = 1 pure super Yang-Mills (SYM) theories in three [6J 
and four dimensions [7], some lattice models are shown to be free from any fine tuning 
at least perturbatively. (Lattice models for three- and four- dimensional SYM theories 
are constructed by orbifold or twisting methods preserving some of supersymmetries [5], 
although they require fine tuning in taking the continuum limit.) In order to overcome 
this difficulty, one possible direction is to pursue new discretization methods different 
from conventional lattice. 

For one-dimensional theory (matrix quantum mechanics) a powerful "non-lattice" 
technique [8] is applicable. Maximally supersymmetric matrix quantum mechanics has 
been studied extensively, and remarkable quantitative agreement with the gauge/gravity 
duality conjecture has been obtained [9| [TP]. (Qualitatively consistent results are ob- 
tained also from lattice simulation [H].) □ For two-dimensional J\f = (2,2) SYM, non- 
perturbative evidences for the lattice model presented in [H] to require no fine tuning 
have been given by numerical simulation for the gauge group G = SU (2) in [15] and for 
G = SU(N) with iV = 2,3,4,5 in [16] . Furthermore, [TT] has shown that the model 
constructed in [18] is free from the sign problem and gives the same physics as that in 
[H] after an appropriate treatment of the overall U(l) modes. M = (8, 8) theory has also 
been simulated in [19] in order to study the black hole/black string transition. (For other 
numerical studies in the context of the gauge/gravity duality, see e.g. [201 EI]-) Com- 

4 In particular the simulation results are consistent with the existence of the threshold bound state [TU] , 
which is an important ingredient of the Matrix theory conjecture [1] . Theories with less supersymmetry 
have also been studied [T2] and the result strongly suggests the threshold bound state does not exist in 
those cases as expected from the calculation of the Witten index. For SU(2) theory with four supercharges 
the energy spectrum has been studied in |13j by using the Hamiltonian approach. The simulation results 
[12] look consistent with the spectrum calculated in [13] . 
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bining the non-lattice or lattice techniques with matrix model/non-commutative space 
techniques [22| 123] . three-dimensional theory can be obtained as a theory on fuzzy sphere 
[24] . Also, in the planar limit, four- dimensional theory can be obtained using a novel 
large- N reduction technique [25j|26] inspired by the Eguchi-Kawai equivalence [27]. How- 
ever, four- dimensional theories of extended supersymmetry at a finite rank of a gauge 
group were out of reachjf] 

Recently we proposed a new regularization scheme for four- dimensional M = 4 SYM 
with G = U(k), which is free from fine tuning [2H1 EHj- It is a hybrid of two-dimensional 
lattice [311 HU [32] (see also [331 CES EU EHJ ESI E0 EE]) and matrix model techniques [21]. 
We regularized a two-dimensional SYM with plane wave deformation, which has a fuzzy 
sphere classical solution, and two additional non- commutative dimensions are generated 
by the fuzzy sphered In four-dimensional N = 4 theory, the commutative limit of non- 
commutative space is believed to be smooth [HJ H2J 03] . Therefore it is expected that we 
can numerically simulate four-dimensional M = 4 SYM on R 4 using this formulation. 

In this paper we provide a non-perturbative formulation of two-dimensional M = (4, 4) 
SYM and four- dimensional M = 2 SYM on non-commutative space, which is analogous 
to the one given in [29] . We first deform the action of two-dimensional M = (4, 4) 
SYM by a mass parameter M with keeping all the supersymmetry. The gauge group is 
U(N) or SU(N). As a result of this deformation, flat directions of three scalar fields 
are stabilized, and not only the trivial configuration but also fuzzy sphere configurations 
become supersymmetric classical solutions. We further introduce another mass parameter 
m, which keeps two supercharges, Q±, in order to lift up the flat direction of the remaining 
scalar field. Then the deformed theory is formulated on a two-dimensional square lattice 
in a manner to preserve Q± exactly. Here we use a prescription developed by one of 
the authors, F.S. [31] [II]- Thanks to the deformations by M and m, we can solve the 
problem of the running of the vacuum expectation values of the scalar fields along the flat 
directions. In this sense, the formulation in this paper can be regarded as a modification 
of the formulation in [T3] to stabilize all the flat directions of scalar fields with keeping Q± 
supersymmetries. We will give a perturbative argument that the continuum limit of the 
two-dimensional lattice does not require any fine tuning. Note that the deformation by 
m does not affect the fuzzy sphere solutions. Namely, they are still solutions preserving 
Q± in the lattice theory even after introducing m. We next consider the lattice theory 
with the gauge group U(N) expanded around a ^-coincident fuzzy sphere background. 

5 Number of fine tunings in the lattice model has been estimated in [28] . 

6 Refs. ED] discuss similar construction of four-dimensional non-commutative spaces from zero- 
dimensional matrix models. 
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Taking the lattice continuum limit first, we obtain four-dimensional M = 2 SYM on 
R 2 x (Fuzzy S 2 ) with the gauge group U(k) deformed by the parameter m. Here the matrix 
size is given by N = kn, the radius of the fuzzy spheres R — -A, the noncommutativity 
parameter 9 = jjs-, the four- dimensional gauge coupling g\ d = 27r0g| d , and the naturally 
introduced UV cutoff ^(n — 1). Although the preserved supercharges at this stage are 
only Q±, the supersymmetry breaking is soft by the mass parameter m. Therefore all 
the supersymmetry is recovered by simply taking the limit of m — > 0. Finally, we take a 
large- N limit sending M — > with fixing k and 0. In this limit, the fuzzy sphere becomes 
non- commutative Moyal plane R@. Because of the full supersymmetry preserved upon 
turning off M, we can strongly expect that four- dimensional M = 2 SYM on R 2 x R| is 
realized without any fine tuning. 

This paper is organized as follows. In the next section, we review continuum M = (4, 4) 
SYM theory in two dimensions and rewrite the action in the so-called balanced topological 
field theory (BTFT) form. In the section 3, we add appropriate terms depending on the 
parameter M to the action, so that all the supercharges are preserved. Furthermore, 
the mass m is introduced to stabilize all the flat directions. In the section 4, we put the 
deformed theory on a two-dimensional lattice in a manner to keep two supercharges. In the 
section 5, we present an intriguing scenario leading to four-dimensional M = 2 SYM with a 
finite-rank gauge group U(k), starting with the two-dimensional lattice formulation given 
in the section 4. The section 5 is devoted to conclusion and discussion. In the appendix 
A, we explain how to derive the mass deformation by M that keeps all the supercharges. 
In the appendix B, the explicit form of the lattice action is presented. 



2 Continuum two-dimensional Af = (4, 4) supersym- 
metric Yang-Mills theory 

In this section we recast the two-dimensional M = (4, 4) SYM theory on (Euclidean) R 2 
with gauge group G = U(N) or SU(N) to a convenient form for our lattice formulation 
(BTFT form [B]). The action of the theory reads 

S2d= wJ d2xTl t^ 2 + \ {v " xl)2 ~ i [x/ ' xJ]2 

+ ^ T (V X + l2 V 2 )V + ^ T 7 /[X 7 , *]] , (2.1) 

where /i = 1, 2, / = 3, • • • , 6, = <9 M + i[A^, •], and all the fields are in the adjoint repre- 
sentation of G. They are expanded by a basis of the representation T a (a = 1, • • ■ , dim(G)) 



3 



as = A a ^T a , ■ ■ ■ . The fermion \1> is an 8-component spinor which is real, \l/ a * = \I/ a , 
but not Majorana. The gamma matrices % (i = 2, ■ ■ ■ ,Q) are 8x8 matrices satisfying 
{7i, 7j} = — 25ij and 72 • • • 76 = —Us- Their explicit form we use is 
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(2.2) 



where 0x2,3 are Pauli matrices. This theory preserves eight supercharges and the super- 
symmetry transformation of the fields is given by 



8'A 1 = e T #, S'A 2 = e T 7 2 #, 5'X 1 = e T 7/ *, 



ST' Vl 



6'V 



-F1272 - (^iX 7 ) 7/ + (V 2 X T ) l2I + -[X / ,X J ] 7/J ) e, 



(2.3) 



where e is an 8-component supersymmetry transformation parameter. 

As a preparation to construct a lattice theory later, we transcribe the theory in terms 
of topologically twisted variables. The two-dimensional M = (4, 4) SYM has three global 
U(l) symmetries: the rotation of (x l , x 2 )-plane U(1)e, the R-symmetry U(1)r, and an- 
other rotation U(l)v whose origin is the chiral rotation in four-dimensional M = 2 SYM 
from the viewpoint of dimensional reduction. We write the fermion \1/ as 



^ ~~ ^2 (£ R , ( R , Ci, £r, £li Cri Cl) 5 
and define the following complex combinations, 

A fl =&+<& x R =e R -ici x L =e L +icL \ L =e L -Ki 



(2.4) 



(2.5) 



The [/(l)-charges of the fields are summarized as 
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where U(l)d is the diagonal subgroup of U(1)e x U(1)v- We rename the fields based on 
the symmetry U(l) d as 

B = X 3 , C = 2X\ (f)± = X 5 ± iX 6 , 
As = ^ + #-2) , A L = -J= ^-x- + 



1 1 / i 

^R = ^= (V'+l + #+2) , l/> L =^= ( X+ + 



Correspondingly we define a new expression of the 8-component spinor, 

^ (0) = (V+i, ^+2, x+, 4>-i,i>-2, X-, ^v-^j , 
which is related with ^ by a unitary transformation, 

= u s ¥°\ 



(2.6) 



(2.7) 



(2.8) 
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(2.9) 



Next we consider the two supercharges, Q' + and Q'_, corresponding to the supersym- 
metry parameters e + and e_, 

(2.10) 



with 



e + = 
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(e± : Grassmann numbers) 



respectively. The transformation of the fields by Q'± is given by 

~ 1 i 

Q'± H n = [4>±,tl>*A =F 2^'^] =F 2 P ^±. 

Q±5 = x±, Q' ± X± = ±[0 ± ,5], Q' tX ± = \[C,B]tH 1 

Q'±H = [0±,Xt] =•= 7^[ B iV±] =F ^[C,X±], 

Q' ± C = 77±, Q±?7 ± = ±[0±, C], Q' T V± = T[<f>+, 0-], 
Q±0± = 0, Q' T 0± = T?7±, 

where H and if M are auxiliary fields. It is seen that Q'± satisfy 

Q'l = (infinitesimal gauge transformation with parameter + ) , 



(2.11) 



(2.12) 
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Q' 2 = (infinitesimal gauge transformation with parameter —4>-) , 
\Q'+,Q'-} = (infinitesimal gauge transformation with parameter C) . (2-13) 

Then we can rewrite the action (12. ip in the so-called BTFT form [H]; 

S 2d = Q'+Q'-F, (2.14) 



with 



T = / d z xTr 
did 



-iB<& - ip+f.tp-a - x+X- - 



(2.15) 



where $ = 2Fyi- After integrating out the auxiliary fields, the explicit form of the action 
reads 



1 

did 



5 2 d =^r I d 2 xTr 



F 2 2 + (V^B) 2 + \{V,Cf + V^+Vrf 



+ \[<P+, <P-] 2 + [B, B] - Cf + I[<7, 0+] [0_, C] 

+ 2i X - (PiV>+2 - V 2 ^ +1 ) - 2i X+ - W-i) 

+ 2B({V+i,'0_ 2 }-{V+2,'0- 1 }) 

- + X+[0->X+] +x-[c,x+] - x+t- 8 ^-] - X-t-B^+l 



(2.16) 



This action is manifestly symmetric under an SU(2)r subgroup of the R-symmetry 
group 577(4). The generators of the SU(2)r are represented as 



-74 
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+ 



5 



SC a (x) 
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which satisfy the SU(2) algebra, 

[Jo, J±±] = ±2J±±, [J++, J — ] = Jo- (2.18) 

We see that (4'+^, (v+i ~ V-) an d (Q' + ,Q-) transform as doublets and 

((/)+, C, —0-) as a triplet under the SU(2)r transformation. 



3 Mass deformation with keeping 8 supercharges 

We next deform the action ( 12.1 ft by introducing a mass parameter M; 

S 2 d,M = S 2 d + Sm, (3.1) 



with 



aid 



'aid 



d xTi 

d 2 xTr 
2M 



1 (M 

2 \6 
M 2 fl 

T 



p\2 



M 



C(x) 2 + <P + (x)<P_(x) - —C{x) [4> + {x), 4>-{x)\ 



, . . „ , . 2M , . . , . . 2M . . . . M / \ / \ 
+ i—B(x)F 12 (x) + —ip+^xjip-^x) + —x+{x)x-{x) - —ri + (x)ri-{x) 
6 6 6 b 

(3.2) 

where p,q,r — 4, 5, 6. This deformation is derived from an eight-supersymmetry analogue 
[45] of the plane wave matrix model [22]. The derivation is summarized in the appendix 
lAl As discussed there, ( 13 .ip is invariant under the supersymmetry transformation, 



5 = 5' + 5 



M 



with 5' given by (12. 3p and 



5 M A^ = SmX 1 = 0, 5m^ 



M 



A p 7 P 7456e. 



(3.3) 



(3.4) 



Namely, the deformed theory (13. 1 p still preserves eight supercharges. 

In order to rewrite (13. ip in the BTFT form, we define the deformed supercharges 
Q± through the deformed supersymmetry transformation (13.31) and the supersymmetry 
parameters ( 12 . 1 01) . The Q± transformation of the fields is E 

Q±A IX = ip ±fl , Q ± tp ±fl = ±W fl <f)±, <2 T V±,u = ippP =f 



7 The transformation of the auxiliary fields is determined so that relations (|3.6|) hold. 
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Q±5 = X±, Q±x± = ±[0±,5], g TX± = i[C, J B] T i/, 

= [0±,X T ] ± ^[S,77±] =F i[C, X ±] + yX±, 

2M M 
<5±C = ?7±, Q±^± = ±[0±, C] + — 0±, Q T r/± = =f[0+, ± yC, 

Q±0± = 0, Q T (P ± = TV±- (3-5) 
We can check that Q± satisfy the nilpotency relations, 

Q 2 + = (infinitesimal gauge transformation with parameter + ) + — J++, 

Q 2 _ = (infinitesimal gauge transformation with parameter — 0_) — —J , 

M 

{Q + ,Q-} = (infinitesimal gauge transformation with parameter C) — —Jo- (3.6) 



Using Q±, the action (13.11) is expressed as @ 



S 2d ,M = (Q+Q- - yj T, (3.7) 

where T is identical with ( 12. 15ft . Although Sid,M is not precisely Q+Q--exact, it is Q±- 
invariant. In fact, since T is gauge and SU(2)r invariant, 

J±±F = JoJ 7 = 0, (3.8) 
and (Q + ,Q-) is a doublet of SU(2)r, 

J±±Q t = Q±, JoQ± = ±Q±, (3.9) 



we sec 



Q+S 2 d,M — Q\Q-J~ — -^-Q+J~ 

Q-S 2d>M = ({Q+, Q-}Q- - Q+Q 2 -) J" - ^Q-F 

MM M 
= -y JoQ-F + y Q+ J— J 7 - —Q-J 7 = 0. (3.10) 



This kind of deformation is discussed for various SYM models in 



In the next section, we consider a two-dimensional lattice theory corresponding to 
this theory. If this theory is naively put on a lattice, however, we soon find that it is 
hard to perform a numerical simulation because of the flat direction along B(x), which 
causes running of the scalar field. In order to avoid it, we further deform the theory by 
introducing an additional mass term to T\ 



T ->• T + A.F, 



(3.11) 



with 



AT 



aid. 



d 2 xTr 



rn 



B(x) 



(3.12) 



where m is a real constant. This deformation clearly preserves the supercharges Q±. After 
these deformations by M and m, the action becomes 



iM,m 
5 2d 
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Q+Q- - y ) (? + &n 



(3.13) 



'aid . 



d xTr 



F 2 



1\ Bf-\{Vt 



\[<P + A-f + [BA + \[^B\ - \[B,C] 2 + i[C,0 + ][0_,C] 
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+ — ( -c( x y + m^-(^) ) - —c{x) [4 + (x), <p.(x 



M 



m ( M m 
-77 Hr + IT B « 



+ 2i X - (V^ +2 - V 2 ip +1 ) - 2i X+ (Pi^-2 - W-i) 
+ 2S({^ + i,V- 2 }- 

- 2^_^_^0+ - 2^+^+^- 
111 



We see that m must satisfy 



2M 



< m < 0, 



(3.14) 



(3.15) 



in order for B(x) to have a positive mass squared. 
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Looking at the bosonic part of this action, we see that there are two types of classical 
solutions; the trivial solution, 



C(x) = (f>±{x) = B{x) = 0, (3.16) 



and the fuzzy sphere solution, 



2M M 
C{x) = —L 3 , ± (x) = —{In ± iL 2 ), B(x) = 0, (3.17) 

where L a (a = 1, 2, 3) belong to an iV-dimensional (not necessary irreducible) representa- 
tion of 577(2) generators satisfying [L a , L 6 ] = ie a b c I c - Around these solutions, there is no 
flat direction because of the mass terms. Also, the Q± transformation of rj±, 

2M M 
Q±V± = ±[0±, C] + — ± , Q tV ± = T[0+, 0-] ± yC, (3.18) 

shows that these solutions preserve the Q± supersymmetry. We here emphasize that the 
shift of m does not affect the fuzzy sphere solution. Furthermore, this solution preserves 
all the eight supercharges of the continuous theory in the limit of m — > 0, as seen from 
the supersymmetry transformation (1A.29j) at m = 0. 



4 Lattice formulation for two-dimensional J\f = (4, 4) 
supersymmetric Yang-Mills theory 

In this section we put the deformed theory on a two-dimensional square lattice with lattice 
spacing a. In this formulation, the supercharges Q± are preserved on the lattice, the gauge 
field is expressed as a link variable U^(x) = e 1 ^^ e G as usual lattice gauge theory, and 
all the other lattice fields are defined on sites and are made dimensionless by multiplying 
suitable powers of a to the continuum counterparts: 

(scalars) lat = a(scalars) cont , (fermions) lat = a 3/2 (fermions) cont , 

(auxiliary fields) lat = a 2 (auxiliary fields) cont , = a 1/2 Q c £ nt . (4.1) 

Also, dimensionless coupling constants on the lattice are 

9o = a92d, M = aM, m = am. (4.2) 

The supersymmetry transformation is realized as 

Q±Up(x) =i^±^(x)U^(x), 
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Q±Hf,{x) = - ~ [ip T ^(x), (f)±(x) + Ufj,(x)(f)±(x + fx)U^(x) ] ] 

± ~ frM*), C(x) + U,(x)C(x + fi)U^) ] ] 

i 1 
=F 2 V nV±(x) ± -[ifj ±t ,(x)tfj ±t ,{x),tp Tfl {x)} 



% 

+ 2 



ip±^x),H^(x) +^?p±n(x), (4.3) 



for the lattice fields U^x), ip±ii{ x ) an d H^x), and transformation of the other fields is the 
same as the one in the continuum theory ( 13.51) with the obvious replacement M — > Mq. 
Here we have used as a covariant forward difference operator, 

2V4(x) = U^x)A(x + flU^x? - A(x), (4.4) 

for any adjoint field A(x). In order to construct a corresponding lattice action, we take 
the lattice counterpart of <3> as 



•w = W*> s i-4ii-K,wir (4 - 5) 



for G = C/(JV) and 



<S(x) = $5C/(JV)(a;) = $[/(jv)0) - ^Tr ($ u{N) (x)^ 1 N , (4.6) 

for G = SU(N). Here U flv (x) = U fM (x)U u (x + fi)U^(x + z>)^ I/ (x)^ is a plaquette variable, 
e is a constant satisfying 0<e<2forG = £/ (N), and < e < 2 a/2 for iV = 2, 3, 4 and 
< e < 2-y/iV sin(7r/iV) for N > 5 for G = SU(N), and the norm of a matrix is defined 
by \\A\\ = v/Tr(A4t) [H]. 

We then put the two-dimensional theory (I3.13P on a lattice using the same form of 
T and AJ 7 in (12.151) and (13.121) together with the trivial replacement -4- / d 2 x — >■ -j 
M — )■ Mo and m — >■ mo- The obtained Q-t-invariant lattice action is 

■ (Q+Q- ~ it) (^iat + AJi at ) , ||l-[/ 12 (z)ll < e for v a; 

>->lat = \ (4.7) 

oo, otherwise 



with 



■Flat = -o 5^ Tr 
00 ~ 



-iB(x)$(x) - ip + ^{x)i)_^(x) - x+(x)x-{x) - ^7]+(x)r]-(x) 
12 



The explicit expression of the lattice action is given in the appendix El 



(4- 



4.1 Absence of flat direction and realization of the physical vac- 
uum 

Let us check that the lattice action has the minimum only at the pure gauge configuration 
^12(20 = ljv which guarantees that the weak field expansion UJx) 



1 + iaApix) 



+ 



A 

2! A 4 



» 2 + 



is allowed in the continuum limit so that the lattice theory converges 
to the desired continuum theory at the classical level. After integrating out the auxiliary 
fields, bosonic part of the action S\ at takes the form, 



5, 



lat 



So 4" 



-T (f + ! Y ]B(X) 



+ i 



Mn mo . , 

+ Ui.r)<\H.r) 



PDT, 



(4.9) 



where Spot denotes positive (semi-) definite terms given by HB.5[) . We will treat the second 
term, which is purely imaginary, as an operator in the reweighting method, and consider 
the minimum of the remaining part of . If the condition — ^y 2 - < rriQ < 0, the lattice 
counterpart of (I3.15p . is satisfied, the mass terms in (14. 9 p fix the minimum at 



B(x) = 0, 

which is independent of Spot- At this minimum, Spdt becomes 
Spdt 4E Tr fe^V^)) 2 + (i[X*(x),X«(x)] + 

-(U 12 (x) - U 21 (x)) 2 } 



(4.10) 



X, 



2n 



4 * 2 4- (1 - ^ 



(4.11) 



\1-Un(x)\\ 2 ) 
with ( 12. 6 p for p,q,r — 4, 5, 6. 

Looking at the first line, we see that the trivial solution ( 13 .161) and the fuzzy sphere 
solution ( I3.17P are still classical solutions of the lattice theory by taking into account the 
replacement M — > Mq. In the same manner as in the continuum theory, there is no flat 
direction around the solutions; we can perform a stable numerical simulation with keeping 
two supercharges. Note that the fuzzy sphere solution plays a crucial role to discretize 
four-dimensional Af = 2 SYM in the next section. 

As discussed in [TJ], in the last term of ( 14. lip representing the gauge kinetic term, the 
admissibility condition singles out the trivial minimum Ui2(x) = In- It shows that the 
lattice action has a stable physical vacuum and unphysical degeneracies of vacua do not 
appear. 
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4.2 Absence of the species doubler 



Let us confirm that there is no species doubler in the kinetic terms of this lattice action. 
Setting Upix) = 1, the kinetic terms for bosons and fermions become 

S ( 2 B) =\ E Tr {(VM*)) (A^-0*)) + 7 (^C(x)) 2 + {A,B{x)f 
9o x I 4 

+ f + icW 2 ) - f (f + f ) (4.12) 

S< 2 F ' -^E^ (* l0,TG ^ ( a « + a ;) * !0) + * ,0,T -P^ (a, - a;,) *i o > + *(»> T x«i»(">), 

(4.13) 

respectively. Here A M and A* are forward and backward difference operators; 

A,f(x) = f(x +p)- f(x), A;fix) = fix) - fix - /}), 



(4.14) 



the matrices and P^ are given by 

— «72 



icy 2 



( 



G 2 = i 



Pi = i 



J 

ia 2 \ 



\-0"3 



-CTi 



0"! 



\io-2 

and the mass matrix .M is 



P 2 = * 



1 2 \ 



0-3 



-cr 3 



(4.15) 



-m d 



m d \ [M M M m M 

^ = d ia g(—,— , — + -,- — 



(4.16) 



Note that G M and P M are anti-hermitian matrices and hermitian matrices, respectively, 
satisfying 

G v } = -28^, {P„ P v } = 25^, {G„ P v } = 0. (4.17) 

The bosonic part ( 14. 12ft takes the form of the standard lattice kinetic terms of bosons; 
no doubler appears in the bosonic sector. For the fermionic part (14. 131) . the kernel in the 
momentum space takes the form, 

2 



T>Fip) = £ [^^W - 2P ^in 2 (^) 



(4.18) 



14 



at Mq = mo = 0. Since the mass terms have the same structure as in the continuum, it is 
sufficient to consider the kinetic terms without the mass terms for our aim. Using (I4.17p . 
we can easily see 

V F {pf = j^Asm 2 . (4.19) 

Since T>p(p) is hermitian, it shows that only the origin (p±,P2) = (0,0) gives the zero of 
T> F (p) in — - < Pn < that is, there is no species doubler in the fermionic sector as well. 

4.3 Absence of parameter fine tunings 

Next, we discuss in the perturbation theory that the desired quantum continuum theory 
is obtained without any fine tuning. In the theory near the continuum limit with the 
auxiliary fields integrated out, let us consider local operators of the type: 

O p (x) = M™ V {x) a d^{x) 2 \ p = m + a + P + 3-f (4.20) 

where (p(x), ip(x) and d denote bosonic fields, fermionic fields and derivatives, respectively. 
M represents M or m. The mass dimension of O p is p and m, a, (3, 7 = 0, 1, 2, - • • . 

From dimensional analysis, radiative corrections from ultraviolet (UV) region of loop 
momenta to O p have the form, 

(-Jr c o a p ~ 4 + ci a?' 2 + g 2 2d c 2 a p + ■ ■ • ^ j d 2 xO p (x), (4.21) 

up to possible powers of ln(aM). Co, ci,cs are dimensionless numerical constants. The 
first, second and third terms in the parenthesis are contributions from tree, 1-loop and 
2-loop effects, respectively. The "• • • " is an effect from higher loops, which are irrelevant 
for the analysis. 

Since relevant or marginal operators generated by loop effects possibly appear from 
nonpositive powers of a in the second and third terms in f)4.2ip . we should look at operators 
with p = 0,1, 2. They are tp, Mip and ip 2 except for non-dynamical operators like 1, M, M 2 
and d<p. For G = U (N), although only the candidates for ip is Tr B from gauge and SU{2)r 
symmetries, it is not invariant under Q± supersymmetries; it is forbidden to appear. 
Similarly, Mtp and ip 2 are not allowed to be radiatively generated. For G = SU(N), we 
may consider ip 2 alone, whose candidates are not generated by the symmetries. 

Therefore, in the perturbative argument, we can conclude that any relevant or marginal 
operators except non-dynamical operators do not appear radiatively, meaning that no fine 
tuning is required to take the continuum limit. In particular, if we consider the lattice 
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theory around the trivial minimum C = <p± = 0, the mass-deformed two-dimensional 
M = (4, 4) SYM is obtained without any fine-tuning. Also, after taking the limit of 
m — > 0, we can safely take the limit of M — > to reach the undeformed theory because 
of the exactly preserved eight supersymmetries. Thus we can use this lattice theory as a 
non-perturbative definition of two-dimensional M = (4, 4) SYM theory. 

5 Four- dimensional J\f = 2 super symmetric Yang-Mills 
theory in the non-commutative space 

In this section, we discuss a scenario to obtain four-dimensional M = 2 SYM on M 2 x 
(Fuzzy M 2 ) from the lattice formulation given in the previous section. 

Let us consider the lattice theory for G = U (N) around the minimum of fc-coincident 
fuzzy sphere solution (I3.17P with M replaced by M , 

L a = <g> l fc (a = 1, 2, 3) and N = nk. (5.1) 

are generators of an n(= 2j + l)-dimensional irreducible representation of su(2) 
corresponding to spin j. 

First, we take the continuum limit of the two-dimensional lattice theory. Then, we 
obtain four- dimensional M = 2 U(k) SYM on R 2 x (Fuzzy S* 2 ) deformed by the mass 
parameter m. The fuzzy S 2 has the radius R = -Jj and its noncommutativity is char- 
acterized by the parameter = Tp-- The UV cutoff A is naturally introduced by the 
size of the matrix; A = 4p • 2j. Although these properties of the fuzzy S 2 can be seen by 
doing a similar calculation as presented in Refs. [47J [2U [25], let us give a brief argument. 
Momentum modes of a field, say B, on two dimensions are expanded further by fuzzy 
spherical harmonics: 

2j J 

%) = E E (5-2) 

J=0 m=-J 

corresponding to the expression (15.11) . The fuzzy spherical harmonic Yj*m is an n x n 
matrix whose elements are given by Clebsch- Gordon coefficients as [25] 

W = v^ E (-V^'Cjrl-r* \jr)(jr'\ (5.3) 

r,r'=-j 

with an orthonormal basis \j r) representing La i n the standard way: 
(4 n) ±i4 n) ) \jr) = y/(j*r)(j±r+l) \jr±l), 
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LT\jr) = r\jr), (5.4) 

and the modes bj m (q) are k x k matrices. It is seen that the fuzzy spherical harmonics 
are eigen- modes of the Laplacian on the fuzzy S 2 : 

t (f ) V?u4 B \*j»] = (f ) 2 j{j+i)m (s.s) 

giving the rotational energy with the angular momentum J on the sphere of the radius 
R = jj. The UV cutoff A = ^ • 2j can be read off from the upper limit of the sum 
of J in the expansion (I5.2p . The fuzzy S 2 is a two-dimensional non-commutative space, 
which is analogous to the phase space of some one-dimensional quantum system, and the 
noncommutativity to the Planck constant h. The quantum phase space is divided into 
small cells of the size 27rh, whose number is equal to the dimension of the Hilbert space. 
Correspondingly, the area of the S 2 is divided into n cells of the size 2tcQ: 

4nR 2 = n ■ 2tt0, (5.6) 

leading to the value = -M-- 

As stressed in the previous section, the supersymmetry is softly broken from eight to 
two because of the mass parameter m at this stage. The eight supercharges are recovered 
by taking the limit of m — > with fixing M. 

Next we take the limit of n — > oo with fixing and k. In this limit, M and A become 

M oc n~ 1/2 ->• 0, A oc n 1/2 ->• oo, (5.7) 

and the fuzzy S 2 is decompactified to the non-commutative Moyal plane M@. Since the 
fuzzy S 2 solution preserves eight supercharges after taking m — > 0, it is strongly expected 
that the theory becomes M = 2 U(k) SYM on M 2 x M|, after taking the above limit. The 
gauge coupling constant of the four- dimensional theory is given in the form 

gl = 2nQg 2 2d . (5.8) 

After taking this limit, the expansion (15. 2p by the fuzzy spherical harmonics can be 
essentially transcribed to the one by plane waves on M@: 

B{q) = J (02 e ^®%)' 

where q and q are two-momenta on M 2 and M@ respectively, the position operator x = 
(21,22) on Rq satisfies [21,22] = i®, and q = (q, q) represents a four-momentum. The 
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modes 6(q) in the four- dimensional space are k x k matrices. It is easy to calculate the 
inner product between plane waves on R@: 

Tr (e^V**) = ^ 5 2 (p + q), (5.10) 

which leads to the 9-dependence of the relation (15. 8p . 

Although more investigation is needed to clarify whether the @ — > limit of the 
theory is continuously connected to the commutative four-dimensional M = 2 SYM on 
M 4 or not§|, to the best of our knowledge this formulation gives the first non-perturbative 
formulation free from fine tuning for four- dimensional SYM with eight supercharges. 

6 Conclusion and discussion 

In this paper, we deformed two-dimensional M = (4, 4) SYM theory with the gauge 
group U(N) or SU(N) by a mass parameter M with preserving all supercharges and 
expressed the deformed action in BTFT form. We further deformed it by introducing 
an additional mass parameter m in a manner to keep two supercharges, Q±, in order to 
lift up all the flat directions of the scalar fields. We then put the deformed theory on a 
two-dimensional lattice with preserving Q± exactly. The problem of the running of the 
vacuum expectation values of the scalar fields is avoided thanks to the deformation by M 
and m. We also gave a perturbative argument that any fine tuning is not needed in taking 
the lattice continuum limit. Thus this lattice theory around the trivial minimum can be 
regarded as a non-perturbative definition of two-dimensional M = (4, 4) SYM theory. To 
perform actual numerical simulation, it should be checked if the imaginary term of the 
bosonic action (14.91) is managed by the reweighting method, which might cause bosonic 
sign problem independent of the fermionic one. We next considered the lattice theory for 
the gauge group U(N) around a ^-coincident fuzzy sphere solution with N = nk. The 
radius of the fuzzy S 2 is and the noncommutativity of the fuzzy sphere is characterized 
by the parameter 9 = iM-- By taking the lattice continuum limit, we obtained four- 
dimensional N — 2 U{k) SYM on M 2 x (Fuzzy S 2 ) deformed by the mass parameter m. 
It was discussed that, by taking the limit of m — > followed by the limit of M — > with 
fixing k and 0, four-dimensional M = 2 SYM on IR 2 x R|, is realized without any fine 
tuning. 

9 It is naively expected that the 8^0 limit would not be continuously connected to the commutative 
theory because of the ultraviolet/infrared (UV/IR) mixing [48]. There is a discussion, however, that 
non-commutative four-dimensional N — 2 U(k) SYM may flow to the ordinary commutative theory in 
the infrared [49] . 
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In contrast to four- dimensional M = 4 SYM, the commutative limit — > of four- 
dimensional M = 2 non-commutative SYM is expected not to be continuously connected 
to the usual commutative M = 2 SYM because of UV/IR mixing. Even if such expectation 
is true and our scenario does not lead to M = 2 SYM on the usual R 4 , notice that non- 
commutative gauge theory itself is an important subject of research in order to clarify 
non-perturbative aspects of gauge theories. For example, when we consider instantons 
of gauge theories, noncommutativity plays a crucial role to resolve the small instanton 
singularity. In our formulation, we can analyze the dynamical aspect of quite wide class 
of observables of four-dimensional M = 2 non-commutative SYM numerically, which will 
give a strong instrument to reveal the non-perturbative structure of supersymmetric gauge 
theories. 

The deformed two-dimensional theory itself is also interesting on its own. In particular, 
since one can introduce mass terms for all scalars keeping two supersymmetries, and 
hence flat directions (along which all scalars commute each other) are all lifted, one can 
perform stable Monte-Carlo simulation, if reweighting for the imaginary term in (14.91) 
works. (Simulations so far utilized supersymmetry breaking mass terms, which make 
the conclusion more or less obscure.) Moreover, the deformation terms consist of mass 
terms and a Myers term, which are quite similar to the so-called fi-deformation [50J EE] 
which is originally introduced in order to regularize the instanton moduli space of four- 
dimensional M = 2 SYM theoryF°l. In the case of the fi-deformation, the integration 
over the instanton moduli space is localized to discrete points, which makes it possible 
to evaluate the instanton partition function analytically using localization formula in 
equivariant cohomology. On the other hand, the deformation introduced in this paper lifts 
flat directions of the scalar fields, which makes it possible to carry out stable Monte-Carlo 
simulation. It is interesting that a mathematically sophisticated technique like equivariant 
cohomology somehow relates to a technique developed for numerical simulation in this 
paper. It may be a sign that such a mathematical method would give a systematic method 
to construct a non-perturbative definition of supersymmetric gauge theories in the future. 
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A Plane wave deformed two-dimensional Af = (4,4) 
supersymmetric Yang-Mills theory 

In this appendix we explain how to construct the plane wave deformed two-dimensional 
M = (4, 4) supersymmetric Yang-Mills theory. 

A.l BMN type matrix model with 8 supercharges 

Let us start with an eight-supersymmetry analogue jl5] of the plane wave matrix model 



S = R / dt Tr 



^(D t X 1 ) 2 + ^> T D^ + y T T l [X\^] + ~[X\X j ) 2 



(A.l) 



where D t = d t + i[A t , ■ ], i = 2,3,4,5,6, a = 4,5,6 and a' = 2,3. T l are 8x8 real 
symmetric matrices corresponding to —iji, which satisfy 

{r\r J } = 26 ij l 8 , (A.2) 

r 23456 = r 2 ---r 6 = -i. (A.3) 

From this model we construct a two-dimensional theory following [53], by using Tay- 
lor's T-duality. In order to lift the theory to two dimensions, we redefine the fields by a 
rotation on (X 2 ,X 3 ) plane with the angle at as 

J) = U °® ' * = e ^ 3at * (A - 4) 



with 

UJt) = I 

— sin (at) cos(at) 



U a (t) = I C ° s(at) Sin(at) I . (A.5) 
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For the other variables, the hatted variables are the same as the unhatted ones. For 
example, D t = d t + i[A t , ■] = d t + i[A t ,-}. Since T 23 is real anti-symmetric, \1/ T transforms 
as ^> T = ^ T e~5 r23a *. The action in terms of the redefined fields is 



R / dtTr 



2R 2X L ' R 1 4 
2\V6i2/ R 2 ) y 1 AR \ n ) 



dt O t Tr (X 2 X 3 ). (A.6) 



+ ^ X 3 D t X 2 - - — (X a ) 2 - t -^e abc X a X b X c 
R 2 1 2 \3RJ V ; 3R 



a 
~R 

If we set a = ±| and discard the surface term, the mass terms of X a ' vanish, and X 2 
appears only in the adjoint form. Then, the Taylor's T-duality may be performed with 
respect to X 2 . As we will see shortly, compatibility of supersymmetry transformation and 
T-duality singles out a = — |. 

A. 1.1 Supersymmetry transformation 

Supersymmetry transformation of this matrix quantum mechanics is given by 

5X i = itf T Fe(t), 
5 A t = Ri^ T e(t), 

5* = ^(D t x^r+ l -[x\x^ 

+ Ax a r a r 456 - J^_x a 'r a 'r 456 ) e (t) (A.7) 

QR \2R ) 

with 

e(t) = e-^ 5 \ , (A.8) 

where eo is an 8-component constant spinor. This is called "dynamical supersymmetry" 
which is simply referred as supersymmetry in the text. For the case of G = U(N), it is 
also invariant under "kinematical supersymmetry", 

5^ = r](t)l N (A.9) 

with 

V (t) = e^ r45 V (A. 10) 



21 



In terms of the redefined fields, the dynamical supersymmetry transformation becomes 
6X* = ii> T T l e(t), 



8A t 
5$ 



RiV T e{t) 
1 



2R K ' 4 L ' 1 6R 



1 //i 



+ a ) ( X 2 r 3 - X 6 Y 



2R\6 

where e is given by e(t) = e^ v2iat e{t). Here we used T 456 = T 23 to obtain 

^a'-pa'-p456 / \^2-p2 , \V3-p3 \ -p23 \A2-p3 \^3-p2 



rr + x^ = rr j - x d p 



(A.ll) 



(A.12) 



For the case that X 2 appears only in the adjoint form in the supersymmetry transforma- 
tion, Taylor's T-duality procedure keeps the supersymmetry. It uniquely fixes the choice 
of a to 

Then, the parameter e(t) becomes t- independent: 



(A.13) 



e(t) = e 



e(t) 



^(r 23 -r 456 )t 



e — e 0) 



(A.14) 



and the dynamical supersymmetry is expressed as 



SX l 
5^ 



RW T e , 

— (DtX^T 1 + -[X\X j \Y ij + JL X a T a T^ Q \ e . 
2R 4 6R ' 



Also, the kinematical supersymmetry becomes 

W = r)(t)l N 

with 



fj{t) = e-^ V (t) = e^ 56+ ^ 2 % = e^V 



(A.15) 

(A.16) 
(A.17) 



The final form of the action (1A.6|) with the surface term dropped is 
S = R I dt Tr 



^(Aif + ^ T D t i> + i> T r[x\ if] + - A \x\ x j f 



-'—J>' I 'I' - J^X 3 D t X 2 - I f A) 2 (x a ) 2 - iJ!-e abc X a X b X c 

3R 



Jf_,TrTr<23v 

3R 



3R 2 ' 



2 V3.R/ 



(A.18) 
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A. 2 Uplift to two dimensions: Taylor's T-duality 

In order to obtain a two-dimensional theory, we "compactify" the X 2 direction to a circle 
of a radius R as 

X 2 ~ X 2 + 2nR. (A.19) 



By using the argument by Taylor [51], the action ( 1A.18I) is lifted to two dimensions, 
S 



rl/R 

R.H I dt J dcxTx 



J_i£ + * ( A x 7 ) 2 - -(D.X 1 ) 2 
2R 2 ta 2R 2 v ; 2 v ; 



+-m T D t m - m> T T 2 D^ + ^ T r 7 [X 7 , + -[X 7 , X 7 ] 2 
R 4 



ifJi) 2 (x a ) 2 -^^r 23 vi>. 

2 V3i?/ V ; 3R 3R 2 



''' X*F tu -i^-t abc X a X h X c 



(A.20) 



where / = 3, 4, 5, 6. (The hats of the fields were omitted.) 
By setting 



t = -^Xi, A t = RA X , cr = x 2 , 



\/R = L, g = (2L) 1/2 , M = ^ 



R 



(A.21) 



and by rescaling the fermion as \1/ — > we obtain 



S = — I d 2 xTi 
9 



1 -F 2 2 + 1 -(D 1 X I ) 2 -\(D 2 Xr 



+U> T (Dx - t 2 d 2 ) q + i^ T r 7 [x 7 , + i[x 7 , x J ] 2 
-1 (HY { xr - £ *w - m x s Fi2 _ JL^tfr 

2 V 3 / 6 3 3 



where 



d x ■ ■ ■ = 



dxx I dx 2 
o 



(A.22) 



(A.23) 



The dynamical supersymmetry is expressed as 



5Ax 
5A 2 
5X 1 



-^ T r 2 e, 

-z^ T r 7 e, 
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5^ = - 1 F i2 r 2 + (DiX^r 7 + (D 2 X I )T 21 

+ l -\X I ,X J }Y IJ + ^XTT 456 ) e, (A.24) 

where e = — eo / \/2 is an 8-component constant spinor. Note that all the dynamical su- 
persymmetries are preserved. It is in sharp contrast with the case of sixteen supercharges 
[53] . where a half of supersymmetries are broken by deformations. 
The kinematical supersymmetry also remains. It is given by 

6* = f}(x 1 )l N (A.25) 

with 

fj{x x ) = ef r23 ^rj , (A.26) 

where rj is constant. 



A. 2.1 Wick rotation 

In order to obtain the Euclidean action, which is going to be put on a lattice, we perform 
the Wick rotation, 

xi -> -ixi, A x -> iA x . (A. 27) 

The Euclidean action is 



S,- 



— I d 2 xTr 



\f 2 12 + \iD.X 1 ) 2 + V (Dx + l2 D 2 ) * 



+U> T 1I [x I ^}-\\x I ,x J f 

+ Uf) 2 (x°T - + if x*f 12 + ;f< v.V.v' 



where /i = 1, 2, and 7/ = iT 1 , which is identical to ( 13. 1ft . 

Then, the dynamical supersymmetry transformation is written as 



(A.28) 



5 Ax = e T ^, 

8A 2 = e T 7 2 ^, 

SX 1 = e T 7/ ^, 

5* = ( -F 1272 - (DxX 1 )^ + (D 2 X T ) l2I + '-[X 1 , X J ) 1U - f A a 7 a7456 ) e, 



(A.29) 
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and the kinematical supersymmetry is given by 

5'V = fj(xi)l N 

with 

fj(xi) = e l ^~ 723x1 770. 



(A.30) 
(A.31) 



B Explicit form of the lattice action 



In this section, we explicitly write down the lattice action (14. 7p in terms of lattice fields. 
We divide the action into the bosonic and the fermionic parts; 



The bosonic part is given by 



where 
A{x) 
A 1 (x) 



5, 



(B) 
lat 



°lat — «->lat "T °lat • 



H{xf - 2iH(x)A(x) + H^xf + 2iH fl (x)A^(x) 



(B.l) 



+ {V^+ix)) (V^x)) + - {V,C{x)f 

- X - [B{x),C(x)} 2 - [B(x),</> + (x)] [B(x),<f>4x)) 



^B(xf + i ^B(x)H*) 
o 3 



1 - i 
■■-${x) + -moB{x), 

1 1 

: 21-i||l-[/ 12 (x)|P 

x -U 12 (x)B(x) -B(x)U n (x) 

+ U 2 {x - 2)- 1 (B{x - 2)U 12 {x - 2) + U 21 {x - 2)B{x - 2)) U 2 {x - 2) 
x Tr (B{x) (U 12 (x) - U 21 (x))] 

' U-in(x) - U^X 



(B.2) 



2e2 (l-^||l-C/i 2 (^)ll 2 )" 
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! Tr (b{x - 2) (U 12 (x - 2) - C/ 21 (x - 2)) 
{1 - - U 12 (x - 2W) 2 
x (u 2 (x - 2)- 1 (U 12 (x - 2) - C/ 21 (x - 2)) U 2 (x - 2)), 



A 2 (xl 



"21-^||l-C/ 12 (x) 



x 



S(x)i/i2(x) + i72i(a;)S(x) 



- c^x - i)- 1 (c/ 12 (x - i)£(x - i) + B(x - i)u 21 (x - i)) c^x - i) 

, TV (C/ 12 (x) - C/ 21 (x)) N 

' ' U 12 (x) - U 21 (x) 



2e2 (l-^||l-f/ 12 (x)|P) 2 
! Tr (B{x - i) ([/ 12 (x - i) - U 21 (x - 1))) 

+ ^ (i-^||i-t/ 12 (^-i)ll 2 ) 2 

x (c^x - i)- 1 (C/i 2 (x - i) - U 21 (x - 1)) C^x - i)). (B.3) 
Note that A M (x) are hermitian. After integrating out the auxiliary fields, becomes 

where Spdt denotes positive (semi-)definite terms: 



(x) 



+ Sp DT , (B.4) 



•S'pdt — Tr 

1 



U>{xf + (2>^ + (x)) (P^-(x)) + i (^C(x)) 2 



4 [B(x),C(x)] 2 - [S(x),0_(x)] 
+ i [0 + (x),<Mx)] 2 - 1 [C(x),0 + (x)] [C(x),0_(x) 



^C(x) [0+(x), 0_(x)] + ^ Qc(x) 2 + + (x)0_(x) ) + I x (x) 2 + l 2 (x) 2 



In order that the field B(x) has positive mass squared, m must satisfy 

2M 



< m < 0. 



(B.5) 
(B.6) 



The fermionic part is given by 



'lat 



=- VTr 
So 2 ^ 



i^ +M (x)P M ?7_(x) + ^- M (x)2V7 + (x) 
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where 



+ X+ (x) [0- (x) , x+ (x)] - x- (x) [(f)+ (x) , x- (x)] + x+ (x) [C(x) , X- {x)\ 

- V+ (x) [B(x) , x- {x)\ - V- (x) [B(x) , x+ {x)\ 

+ \v+(x)[(f>-(x),r] + (x)} - -ri-(x)[<i> + (x),ri-(x)] - ^rj + (x)[C (x) , r]_(x)} 

- i>+»(x)il>+v(x) ((/>-(x) + U^(x)(f)_(x + foU^x)- 1 ^ 
+ ^(x)^(x) (<j>+(x) + U^x)<p + (x + /i)C^(x)- 1 ) 

- 1 {i; + ,(x),^(x)} (C(x) + U,(x)C(x + W^x)- 1 ) 
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^ 3 + ™oJ X+(x)x-(x) - —r] + (x)ri„(x) 



+ i X -(x)(Q + Hx)) -ix + (x)(QMx)) -iB{x)(q + QMx) 
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Q±Hx) 



-1 



X 



1-^||1-M^)II 2 

- (ip ±1 (x) + C/i(x)V ±2 (x + ijcMx)- 1 ) C/i 2 (x) 

- M*) frM*) + U^x^ix + i)^)" 1 ) 
+ (M*) + U 2 (x)4 ±1 (x + 2)U 2 {x)- 1 ) U 21 (x) 

+ M*0 (V± 2 (x) + C/ 2 (x)V±i(x + ^^(x)- 1 ) 

M^-M*) 1 



+ 
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(l-£||l_M*)ll 2 ) 6 



x Tr 



Ma;) - U 2 i(x) )[V 2 ^ ±1 (x) - V^ ±2 (x) 



Q+QMx) 



g+g_ (u 12 (x)-u 21 (x)) 



fermion 
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1-^||1-M^)II 2 
| i (Ui 2 (x) - U 21 (x)) ^ 1^ 



|1-M^)II 2 
1 



2 £ 2 



g+g_ (c/ 12 (x) + c/ 21 (x)) 
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12 1 6 



(l-^||l-f/l 2 (^)|| 2 ) 

x {-^g + (C/ 12 (x) - C/ 21 (x)) Tr [g_ (C/ 12 (x) + C/ 21 (x)) 
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+ iQ. (U 12 (x) - U 21 (x)) Tr \Q + (U 12 (x) + U 21 (x)) 
i(U 12 (x)-U 21 (x)) 1 



3 e 4 



(l - - U 12 (xW) 

Q+(U 12 (x) + U 21 (x)) Tr Q-(U 12 (x) + U 21 (x)) 



x Tr 



(B.9) 



with 



Q+Q- (u 12 (x)-u 21 (x)) 
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fcrmion 



2 [^+i(x),^_i(a;)] - -C/ 1 (x)[^ +2 (rr + l),V-2^+ l)]^)" 1 

- v + i(x)c/ 1 (x)V-2(^ + i)czi(x)- 1 + V-i(x)c/i(x)v +2 (x + i)^) -1 }*^) 

+ C/i 2 (x){-^ [^ +2 (x), V- 2 (x)] - ^f/ 2 (x)[^ +1 (x + 2), + 2)]f/ 2 (x)- 1 

- f/ 2 W^ + i(x + 2)f/ 2 (x)~V- 2 (x) + C/ 2 (x)V-i(x + 2)U 2 {x)-^ +2 {x)^ 

- {-i foMx), ^_ 2 (x)] - ^U 2 (x)[ij> +1 (x + 2), + 2)]f/ 2 (x)- 1 

- V +2 (x)C/ 2 (x)V-i(x + ^^(x)- 1 + 4_ 2 (x)U 2 (x)if> +1 (x + ^(xy^U^x) 

- u 2 i(x){~ [Mx), V-iW] - ^i(x)[^ +2 (x + i),^_ 2 (x + i)]^^)- 1 

- (^_x(x) + c/i(x)V- 2 (x + i)EMaO _1 )MaO 

x (^ +2 (x) + f/ 2 (a;)^ + i(a; + 2)f/ 2 ( a :)- 1 ) 

+ (^ + i(x) + U^xty+tix + iw^xy^u^x) 

x (^_ 2 (x) + f/ 2 ( a :)^ 1 ( a : + 2)f/ 2 ( a; )- 1 ) 
+ (^- 2 (x) + U 2 {x)^{x + ^^(x)" 1 )^!^) 

x (^ +1 (x) + ^i(x)^ +2 (a: + 

- (iP +2 (x) + U 2 (x)i/j +1 (x + 2)U 2 {x)- l )U 21 {x) 

x (^ 1 (x) + U 1 {x)^ 2 {x+i)U 1 {x)- 1 ), (B.10) 



Tr 



g+g_ (c/i 2 (x) + c/ 21 (x)) 
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+ i (P 2 ^+i(a:)) V-i(^) + \u 2 {x)^{x + ^[/.(x)- 1 (P 2 ^+i(x)) 
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- X - {V^_ x {x)) t/j +1 (x) - hj 2 {x)^ +1 {x + 2)U 2 {x)~ l (P 2 ^-i(a:)) 
+ l -i>+ 2 {x) (V^_ 2 (x)) + 1 (Pi^_ 2 (x)) £/iOr)V +2 (z + l)^!^)" 1 

- ^- 2 (x) (Pi^ +2 (a;)) - ^ (Pi^ +2 (a:)) C/i(x)^_ 2 (x + i)^)" 1 
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(M*) - *7 2 i(aO) Tr |Q_ (C/ 12 (x) + U 21 (x)) 
+ iQ_ (U 12 (x) - U 21 (x)) Tr [Q + (C/ 12 (x) + C/ 21 (x)) 

= | (^+i(x) + c/i(x)v +2 (x + i)^!^)- 1 )^^) 

- C/i 2 (x)(^ +2 (x) + C/ 2 (x)V+i(x + ^[/.(x)- 1 ) 
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x Tr [i (C/ 12 (x) - U 21 (x)) (-P 2 ^_x(x) + Pi^_ 2 (a;)) 

(V>_i(x) + c/i(x)V- 2 (x + i)^)- 1 )^) 
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+ c/ 2 i(x)(^_i(x) + u^x^ix + i)^)" 1 ) 

x Tr 
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= Tr 
x Tr 



i (Un(x) - U 21 (x)) [~V 2 ^i{x) + V^ +2 (x) 
% (U 12 (x) - U 2l {x)) [~V 2 il}-i{x) + Pi^- 2 (x))] . (B.13) 
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